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We investigate the electronic specific heat of overdoped BaFe 2 (Asi_a:Pa :)2 single crystals in the supercon¬ 
ducting state using high-resolution nanocalorimetry. From the measurements, we extract the doping depen¬ 
dence of the condensation energy, superconducting gap A, and related microscopic parameters. We find that 
the anomalous scaling of the specific heat jump AC oc T^, found in many iron-based superconductors, in this 
system originates from a Tc-dependent ratio A/ksT^ in combination with a doping-dependent density of states 
A(ep). A clear enhancement is seen in the effective mass m* as the composition approaches the value that has 
been associated with a quantum critical point at optimum doping. However, a simultaneous increase in the su¬ 
perconducting carrier concentration Us maintains the superfluid density, yielding an apparent penetration depth 
A that decreases with increasing Tc without sharp divergence at the quantum critical point. Uemura scaling 
indicates that Tc is governed by the Fermi temperature Tp for this multi-band system. 


Iron-based superconductors provide new grounds for in¬ 
vestigating unconventional pairing mechanisms in high- 
Tc superconductivity, a long-standing mystery in mod¬ 
ern condensed matter physics. Among the large num¬ 
ber of compounds discovered in this family, the isovalently 
doped BaFe 2 (Asi_a;Pa ;)2 (P-122) system [1] has attracted 
widespread attention. An advantage with the P-122 system is 
that crystals on the strongly overdoped side are clean enough 
to display the de Haas-van Alphen effect (dHvA), yielding 
information about the band structure that can be directly com¬ 
pared with calculations [2-5]. In the system, several meth¬ 
ods indicate the presence of a possible quantum critical point 
(QCP) near optimum doping Xcr — 0.3 [6-10]. In particular, 
the effective mass m* of at least one of the electron pock¬ 
ets at the X-point is enhanced when approaching Xcr [3, 7]. 
Recently, it was reported that the London penetration depth 
Al also has a maximum at Xcr [6]. This is unexpected, since 
this corresponds to a minimum in the superfluid density ps 
when Tc has its maximum, in contrast to the Uemura scaling 
Tc oc nsjm* (x ps = of cuprate and other high-Tc super¬ 

conductors [11, 12]. One of the important pathways to address 
this discrepancy as well as to gain further understanding of the 
pairing mechanism is through specific heat measurements as 
thermodynamic properties are closely linked to the nature of 
the superconducting state. It is known that the thermodynam¬ 
ics of the superconducting state of P-122 is anomalous in cer¬ 
tain ways. The absolute specific heat jump at Tc is found to 
follow the so-far unexplained BNC scaling, AC(Tc) cx T^, 
seen for many iron-based superconductors [13, 14]. How¬ 
ever, thermodynamic evidence of the anomalous behavior in 
Al around Xcr has so far been rather ambiguous [7] and it is 
still debated which anomalies are central to the pairing mech¬ 
anism. 

Here we perform high-resolution specific heat measure¬ 
ments of overdoped BaFe 2 (Asi_a;Pa ;)2 single crystals using 
a nanocalorimetry system capable of providing good absolute 
accuracy on very small samples. From the specific heat data 


we extract the doping dependence of the condensation energy 
AF{0), superconducting gap A(0), density of states at the 
Fermi level N{e'p) as well as microscopic parameters, i.e., co¬ 
herence length ^ab(O) and penetration depth Aab(O) (including 
the related parameters and Al). We find that the anomalous 
BNC scaling may be explained by a joint effect of the dop¬ 
ing dependences of A(0)/fcBTc and A(£f)- While a signif¬ 
icant m* enhancement is observed close to optimum doping, 
it does not lead to a drop in the superfluid density due to a 
simultaneous increase in the superconducting carrier concen¬ 
tration. The evolution of Tc scales with the Fermi tempera¬ 
ture Tp, in accordance with the Uemura relation for cuprates 
and several heavy-fermion systems. This strongly suggests 
that overdoped BaFe 2 (Asi_ 2 ,P 2;)2 belongs to this group of un¬ 
conventional superconductors with an electronic origin of the 
coupling mechanism, governed by Tp. 

High-purity BaFe 2 (Asi_a;Pa ;)2 crystals with x = 0.32, 
0.50 and 0.55 were grown using a self-flux method [14, 15]. 
Small crystals with ^150 p.m side were selected and gently 
cleaved to obtain plate-like samples with shiny surfaces. Spe¬ 
cific heat was measured using a differential, auto-adjusting 
nanocalorimeter operating in constant-phase ac mode to yield 
combined high resolution and good absolute accuracy [16, 
17]. The sample temperature oscillation amplitude was set to 
0.5% of the absolute temperature. Device and thermal adhe¬ 
sive grease addenda were measured separately. To determine 
the absolute specific heat, the crystal volumes were estimated 
from scanning electron and optical microscopy. The result¬ 
ing values at maximum C/T, where the phonon contribution 
dominates, were found within 10% agreement with literature 
values on the parent compound BaFe 2 As 2 [18] as well as non¬ 
superconducting Ba(Feo. 88 ^ 110 . 12)2 As 2 crystals [19]. 

Figure 1 shows the electronic specific heat for three sam¬ 
ples of different composition in magnetic fields H || c. Com¬ 
paring the zero-field temperature dependences of panel (a) 
through (c), there seems to be a tendency of going from a fully 
gapped system for x = 0.32 to a system with pronounced gap 
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FIG. 1. (Color online) Temperature dependence of the electronic specific heat as Ce/T for BaFe 2 (Asi_a;P 2;)2 samples in fields from 0 to 5 T 
for (a) X = 0.32, (b) x = 0.50, and (c) x = 0.55. Field steps are 1 T except in (c) where 0.5 T is also shown. Dotted lines illustrate the 
entropy-conserving construction to determine T^. Inset of (a) shows the total specific heat for the x = 0.32 sample (shown as CjT). Upper 
inset of (b) shows the low-temperature specific heat as C/T vs for the same conditions. The dashed line represents the normal state. Lower 
inset of (b) shows a sample mounted on the calorimeter. Inset of (c) shows the upper critical field pioiTc 2 for the x = 0.55 sample. 


anisotropy and possible nodes for x = 0.55. However, the 
magnetic field dependence also indicates a faster-than-linear 
increase of ^{H) for all samples. This supports the conclusion 
of thermal conductivity and penetration depth measurements 
that possible accidental nodes reside on the high Fermi veloc¬ 
ity electron sheets [ 21 ] to which specific heat measurements 
are not as sensitive. Even for x = 0.32 there is also a sig¬ 
nificant residual specific heat coefficient 7 ^ ^ 5mJ/molK^, 
constituting roughly 17% of the normal state Sommerfeld co¬ 
efficient 7 „. Here, 7 r is largely invariant with doping, in con¬ 
trast to Co-122 where an increase in 7 r away from optimal 
doping can be clearly identified [ 22 ]. 

The phonon background of the x = 0.32 and x = 0.50 
crystals were obtained by applying a simple Debye law and 
entropy conservation on the partially field-suppressed states. 
The normal state of the x = 0.50 sample is illustrated as a 
dashed line in the inset of Fig. 1(b). A 5 T magnetic field is 
sufficient to fully suppress the superconducting state for the 
x = 0.55 sample, as indicated in the bottom inset of Fig. 1(c). 
The Debye temperatures of all three samples are within 10% 
of each other. 

From specific heat, a wealth of information can be gained. 
The results are summarized in Table I together with assump¬ 
tions and relations used. Tc, AC/Tc, and the upper critical 
field slope fiQ\dH^ 2 /dT\ are determined using an entropy- 
conserving construction, as shown in Fig. 1. By integrat¬ 
ing AC'(T) = Ce_s — the condensation energy AF(0) 
is obtained. As seen in Table I, both AC/Tc and AF{0) 
vary strongly with Tc, with higher Tc samples also display¬ 
ing higher AC/Tc and AF"(0). There is a clear increase 
of 7 n when approaching optimum doping, indicating an ef¬ 
fective mass enhancement near Xcr- From 7 n the density 
of states N{eF) is found. An average zero-temperature su¬ 
perconducting gap A(0) can then be estimated from AF(0). 


Both N{eF) and A(0) increase with Tc- The obtained A(0)s 
are comparable to, but somewhat lower than those reported 
by angle-resolved photoemission spectroscopy (ARPES) for 
P-122 crystals with similar compositions [23-25]. Attributing 
AF{0) only to a fraction ( 7 n — 7 r)/ 7 n of the samples would 
increase the average A(0) by 10%-20%, largely eliminating 
the difference between our results and ARPES measurements. 
The ratios A(0)/fcBT'c are in line with measurements on Co- 
122 [ 22 ]. 

AF{0) leads us to the thermodynamic critical field iTc(O). 
The upper critical field i 7 c 2 ( 0 ) is estimated from its slope 
close to Tc and the Werthamer-Helfand-Hohenberg (WHH) 
relation [26]. This estimate gives good agreement with high- 
held transport measurements [ 8 ]. Both critical fields reveal 
significant positive correlation with Tc- Having the critical 
fields, the microscopic superconducting parameters, 5 ab( 0 ) 
and Aab(O) can be attained. They both decrease with in¬ 
creasing Tc, but the decrease in Aab(O) is rather minor, as 
seen in Table I. Note that we make a distinction between 
Aab and Al. Al is the zero-temperature London penetra¬ 
tion depth while Aab = ^ahi/T) is defined at all tempera¬ 
tures through Hc{T). Similarly, we obtain the constant ^0 
from Ginzburg-Landau (GL) relations and relate ^ahi/T) to 
Hc 2 {T) through ^oHc 2 iT) = (I)o/2tt^1^{T). We do not make 
any explicit assumptions on the temperature dependence of 
KciT) = Aab(r)/Cab(r). 

In Fig. 2, we illustrate the variation of several supercon¬ 
ducting parameters as a function of Tc- The present measure¬ 
ments are also complemented with an analysis of the results of 
Ref. [14]. For this analysis, we estimate Hc{0) from the mea¬ 
sured properties near Tc by assuming a scaling Kc(0) oc Tc 
that is roughly seen for the present data, as shown in Fig. 2(a). 
Figure 2(b) suggests that iTc(O) is largely proportional to T^, 
leading to an approximate scaling AF"(0) cx T//. This ob- 
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TABLE L Parameters from specific heat using the specified relations 
and assumptions. 


Parameter 

unit 

X = 0.32 

0.50 

0.55 

Tc 

K 

28.4 

18.2 

12.5 

7n 

mJ/molK^ 

28.0 

20.2 

19.2 

II 

1 

mJ/molK^ 

23.3 

14.0 

13.8 

AC/Tc 

mJ/molK^ 

47.6 

29.2 

15.8 

po|dI7c2/dT| 

T/K 

2.13 

0.96 

0.59 

“AF(0) 

J/mol 

5.16 

1.22 

0.39 

VoiTc(O) 

T 

0.468 

0.230 

0.130 

>o77c2(0) 

T 

44.2 

12.8 

5.38 

‘'iV(eF) 

state s/eV 

11.9 

8.6 

8.2 

=A(0) 

meV 

4.2 

2.4 

1.4 

^Cab(O) 

nm 

2.73 

5.08 

7.82 


nm 

3.15 

5.87 

9.03 

'’Aab(O) 

nm 

182 

200 

229 

'Al 

nm 

156 

165 

212 

-'wf 

km/s 

63.8 

68.0 

60.7 

'‘Tf 

K 

602 

492 

373 


^ from integration of AC. 

^ using A_F(0)/Vin = (0)/^. with Vin from [20], using 122 formula 

unit for mol. 

assuming clean-limit WHH, — 0.73Tc \dH^2/^'^\x ■ 

^ Density of states for both spin directions per 122 cell, using 
A^(£f) = 37n/7r2fc|. 
assuming AF = N(ep) A'^ (0)/4. 

f using/ioJ/^2(0) = '?io/ 27 i-?^ 6 ( 0 ). 

® using fj.o\dH^ 2 /dT\j, = 0o/[27r(O.74^o)^7’c] from the clean-limit GL 
relation C-jf(r) I near Te = (0.74$o) (1 - T/Tc ). 

using/ro-ffc(O) = 0o/2v^7rAab(O)Cab(O). 

‘ using the clean-limit GL relation k.c(Tc) = 0.71Al/0.74^o with Kc{Tc) 
obtained from the ratio of ^^d dHc/dT, with 

VmfJ^o |di7c/dT|^ = AC/Tc {Rutgers relation), 
j assuming ^0 = hv-p/7rA{0). 

^ assuming ep = ^bTf = m*v^/2 

servation is in drastic contrast to the expected scaling 
of conventional superconductors, but follows the same trend 
as previously observed in a number of other Fe-based su¬ 
perconductors [27]. It is consistent with the BNC scaling, 
AC oc T(?, provided that Hc(T) has a normal BCS-like tem¬ 
perature dependence, in which case Rutgers relation would 
give AC/T cx cx [H^{0)/T^]^ cx T^. 

The upper critical field slope HQdH^ 2 ldT, shown on the 
right axis of Fig. 2(b), also reveals a dependence, as re¬ 
ported in Ref. [14]. As a result, ^ab(O) varies approximately 
as Tc as seen in Fig. 2(c). The ratio AC/ 7 nTc, describ¬ 
ing the relative specific heat discontinuity at Tc, is shown in 
Fig. 2(d). Interestingly, this ratio is almost proportional to 
rather than being constant, as it would be in the weak-coupling 
BCS theory (AC/ 7 nTc = 1.43). Thus, 7 n increases with Tc 
when approaching Xci, but not as fast as AC/Tc does. The 
ratio of A(0) to Tc also shows an evident positive correlation 
with Tc, see Fig. 2(e). However, this dependence is weaker, 
going from about 2 for the highest Tc and extrapolating to 
about 1 for Tc = 0. 

From 7 n, the carrier effective mass to* is obtained. Band- 


structure calculations for fully doped BaFe 2 P 2 give the bare 
band electronic specific heat 7b = 6.94mJmol~^K“2 [7], 
From 7n/7b we estimate an averaged effective mass m* jm^, 
taking into account the band mass TOb/TOe of each of the five 
Fermi sheets and their respective contribution to the total den¬ 
sity of states (DOS) [5, 7]. The results are shown in Fig. 3(a) 
together with a logarithmic fitting function of the critical ef¬ 
fective mass enhancement as previously applied to dHvA re¬ 
sults [7]. Two well-documented data points are also included 
in the figure, one for the overdoped BaFe 2 (Aso. 37 Po. 63)2 [4] 
and the other for BaFe 2 P 2 [5], which have been confirmed by 
various separate studies [3, 7]. The systematic vertical shift 
between our data and the previous results reflects different 
normalizations involved. From dHvA measurements [4, 5], 
the 7 hole band is revealed to possess a significantly higher 
effective mass. Hence, to* is shifted upwards when all bands 
are taken into account. 

The relation /toAl = m*juse^ allows the evaluation of 
the superconducting carrier density rig when to* and Al are 
both known. In the inset of Fig. 3(a), rig is plotted as a func¬ 
tion of X. Evidently, Us is displaying a possible vanishing 
around x = 0.65, in the vicinity of which superconductiv¬ 
ity is lost [7, 28]. On approaching Xcr, Us experiences a lin¬ 
ear increase without critical signatures. We thus find that rig 
increases with increasing T^. and condensation energy. Our 
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FIG. 2. (Color online) Approximate scaling dependence of supercon¬ 
ducting parameters on Tc. (a) Kc, (b) He and d77c2/dT, (c) ^ab, (d) 
AC/ynTc, (e) A/^bTc. Filled symbols represent measurements of 
this work whereas partly filled symbols are from analysis of the data 
of Ref. [14]. For the latter, the additional assumption that Kc oc Tc is 
made [crosses in (a)], as approximatively found for the present data. 
Dashed lines are linear fits through origin. 
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FIG. 3. (Color online) (a) Carrier effective mass as a function of 
X. Data points (O) at a; = 0.63 and 1.0 are for the /3 band from 
Refs. [4] and [5]. They are shifted upwards (□) to account for the 
effective mass enhancement of all five Fermi sheets (see text). The 
dashed line is the best logarithmic fit to the effective mass enhance¬ 
ment measured by dHvA, y = 1.0 — 0.91 In (a; — 0.30) [7]. The 
solid line is a similar fit to our data and the two renormalized lit¬ 
erature points, y = 1.87 — 0.70 In (a; — 0.30). Inset: Number of 
superconducting carriers per 122 formula as a function of x. (b) 
Uemura plot showing the approximate linear scaling of Tc with Tp. 
Data points other than those reported in this work were reproduced 
from Ref. [6] and [12]. Tb is the Bose-Einstein condensation tem¬ 
perature for ideal 3D boson gases. Inset: Penetration depth Aab (0) 
as a function of x. 


data suggest a significant deviation from the assumption of 
AC/jnTc = 1.43 used by Walmsley efaZ. [7]. We observed 
a clearly identifiable change in AC/^nTc, up to a factor of 3 
- 4 as a function of x, see Fig. 2(d). This implies that an in¬ 
crease of AC/Tc alone close to a:cr does not necessarily lead 
to a proportional increase of m*, as previously suggested in 
Ref. [7]. Such assumptions may instead lead to an overesti¬ 
mation of Al near Xct- In fact, we find a smooth, continuous 
decrease of Aab(O) with decreasing x, as shown in the inset 
of Fig. 3(b). This is in apparent contradiction with the re¬ 
sults of Hashimoto etal. [6], where Al is reported to display 
a peak around Xcr, increasing roughly by a factor 2 to 2.5. A 
possible explanation for this may be that specific heat mainly 


reflects carriers that are strongly gapped (hole band), while di¬ 
rect penetration depth measurements would weigh more heav¬ 
ily on carriers with long mean free paths but not necessarily 
all paired (electron band). This interpretation would be con¬ 
sistent with stronger inter-band coupling near Xcr, where the 
system is more disordered, and weaker inter-band coupling 
but stronger gap anisotropy for the heavily overdoped, but 
cleaner system. 

Besides m*, specific heat provides additional band- 
structure properties. We estimate the averaged Fermi veloc¬ 
ity Up from A(0) and ^ab(O), see Table I. The wpS obtained 
are of the same order of magnitude as those previously de¬ 
duced through ARPES measurements [29] and change only 
little with X. This is largely supported by ARPES data, where 
a variation in up of less than 25% was observed for a similar 
range of x. In Eig. 3(b) we visualize the variation of as 
a function of Tp (obtained from m* and up) in the Uemura 
plot [6, 12]. Strikingly, our data fall seamlessly into the gap 
formed between cuprates and heavy Eermion compounds on 
the plot, showing an approximately linear scaling relation be¬ 
tween Tc and Tp. This strongly indicates that an electronic 
coupling mechanism may be involved in all these apparently 
very different superconducting materials. We note that this 
scaling relation works best as a function of Tp. Plotted as a 
function of A^^ or rig, P-122 would display a horizontal shift 
due to its multiple bands. 

The Uemura plot here presents a rather different picture 
than that of Hashimoto etal. [6]. We find no anomaly in 
the superfluid density close to optimum doping and see a pre¬ 
served scaling of Tc with Tp. Parameters extracted from spe¬ 
cific heat measurements predominantly represent portions of 
the Eermi surfaces with strong pairing. The non-vanishing 
residual electronic specific heat coefficient, also observed for 
instance in Co-122 [22] indicate that patches of the Eermi sur¬ 
face are only very weakly paired or remain unpaired. While 
regular disorder could possibly cause this, it seems unlikely 
for the case of phosphorous-doped 122, where at least the 
strongly overdoped side is in the clean limit. A more likely 
explanation is a decreasing interband coupling with increasing 
doping X, causing both increased gap anisotropy and weak¬ 
ened superconductivity. 

In conclusion, we find that a higher Eermi temperature is 
beneficial rather than detrimental for the electronically medi¬ 
ated superconductivity in phosphorous-doped 122. With in¬ 
creasing overdoping, superconductivity is suppressed through 
a decreasing superconducting carrier density. 
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